METHOD AND APPARATUS FOR ACQUIRING MULTIDIMENSIONAL SPECTRA 
AND IMPROVED UNIDIMENSIONAL SPECTRA WITHIN A SINGLE SCAN 



RELATED APPLICATION 

This application is a continuation-in-part of PCT/US03/21314 filed 07/07/03, the 
contents of which are here incorporated in their entirety. 

BACKGROUND OF THE INVENTION 

Field of the Invention ^ 

The present invention relates to a new method and apparatus for acquiring uni- 
dimensional and multi-dimensional spectra within a single scan, and more particularly, 
to a method and apparatus for acquiring uni- and multi-dimensional nuclear magnetic 
resonance spectra within a single scan. The present invention also relates to a variety of 
methods employing the principles of these methods or techniques. 

State of The Art 

In terms of applicability to unravel the nature of molecular systems, few analytical 
techniques rival the insight furnished by spectroscopy studies based on the acquisition 
of multidimensional spectra. Such spectroscopies include optical, paramagnetic, 
electron, mass and nuclear magnetic resonance (NMR) spectroscopies. Included 
among this latter category is the method known as magnetic resonance imaging (MRI). 
Whether used to discover new pharmaceutical drugs, to characterize new catalysts, to 
investigate the structure and dynamics of proteins, or to carry out a non-invasive 
imaging diagnosis, few analyses are nowadays performed without the aid of at least 
some form of multidimensional experiment, and especially a multidimensional NMR 
experiment. One drawback of these protocols is that, by contrast to one-dimensional 
spectroscopic methods, multidimensional techniques, and particularly multidimensional 
NMR, requires relatively long measurement times associated with the acquisition of 
hundreds or thousands of scans. This places certain kinds of rapidly changing systems 
in Chemistry, such as proteins changing conformation, analytes flowing through a 
chromatography column, unstable biomolecules, outside their realm. Also unsuitable to 
relatively slow multidimensional NMR analyses are the thousands of compounds that 
can currently be made available by a single combinatorial-chemistry assay. Long 
acquisition times also make multidimensional NMR, as presently performed, ill-suited for 
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in vivo analyses, and for clinical measurements in combination with spatial localization 
techniques. 

SUMMARY OF THE INVENTION 

Accordingly, the object of the present invention is to provide a method and 
apparatus for acquiring multidimensional spectra, and especially multidimensional 
nuclear magnetic resonance spectra, within a single scan. This is accomplished by a 
novel method and apparatus that enables treating a sample to acquire multidimensional 
spectra within a single scan comprising the steps of: (1) partitioning, at least notionally, 
a sample into a set of independent subensembles endowed with different resonance 
frequencies; (2) implementing a polychromatic irradiation of the sample whereby the 
various subensembles are selectively manipulated by a time-incremented series of 
excitation or refocusing sequences; (3) follow these various selective manipulation 
processes by a homogeneous sequence capable of generating an observable spectral 
signal from each of the subensembles; (4) simultaneously monitoring the observable 
signals arising from the various subensembles in a resolved fashion; and (5) processing 
the set of signals acquiring in this manner into a complete multidimensional spectral 
data set. The acquisition of the multidimensional spectra can result from the practice of 
any type of spectroscopy. 

In a further development of the present invention, the object of the present 
invention is to provide a method and apparatus for acquiring multidimensional spectra, 
and especially multidimensional nuclear magnetic resonance spectra, within a single 
scan. For the case of NMR, this is accomplished by a novel method and apparatus that 
enables the acquisition of multidimensional NMR spectra within a single continuous 
scan. Provided that an analyte's signal is sufficiently strong, the invention can, in turn, 
shorten the acquisition time of any multidimensional application or experiment by 
several orders of magnitude. The new invention is compatible with the majority of 
multidimensional NMR pulse sequences hitherto proposed, and can be implemented 
using known magnetic resonance hardware. The same single-scan protocol providing 
the spectroscopic information can additionally be exploited to extract information about 
the spatial location of the analyte's originating the signal. Spatially resolved MRI 
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information can thus be extracted from the protocol at no additional experimental costs 
and/or penalties. 

Additional derivations arise from the partitioning of the sample exploited by this 
invention to collect NMR spectra, that relate not only to multidimensional acquisitions 
but also to the ways in which unidimensional NMR data is collected. One of the 
opportunities enabled by this invention is the use of indirect schemes to detect in a 
single scan the NMR spectra of insensitive nuclei / by monitoring instead the NMR 
signal arising from neighboring nuclei S of higher sensitivity. Another one resides in the 
feasibility of collecting high resolution NMR spectra even upon using inhomogeneous or 
unstable magnetic fields, by manipulating the phases used to excite and detect the 
signals arising from each of the excited slices. 

Accordingly, the foregoing objectives are accomplished by a method and 
apparatus for treating a sample to acquire multidimensional magnetic resonance spectra 
within a single scan (or a small number of scans) comprising the steps of: (1) applying a 
magnetic field gradient on the sample so as to endow spins at different positions within 
the sample with different resonance frequencies; (2) applying a train of frequency- 
incremented radiofrequency (RF) pulses in unison with this gradient (or with an 
oscillating version of thereof), so as to endow spins at different positions within the 
sample with incremented values of their evolution times, thus creating an effective 
spatial encoding of the spins' frequencies (3) applying if needed a homogeneous mixing 
pulse sequence at the conclusion of the various spatial encoding processes, capable of 
creating a set of observable spin signals; (4) capturing the signals thus created from the 
sample while decoding the spins' spatial locations using a second set of acquisition 
magnetic field gradient; (5) subjecting the collected data to a suitable rearrangement 
and Fourier analysis procedure so as to retrieve the final spectrum being sought. 

Lying at the core of this new invention is the application of a magnetic field 
gradient, operated in unison with a train of spatially-selective radiofrequency (RF) 
pulses. This enables the invention in a first step to endow spins at different positions 
within the sample with incremented values of their indirect evolution times, leading to an 
encoding of the internal evolution frequencies of spin coherences along a spatial 
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coordinate. This indirect-domain information is preserved throughout a mixing process 
step, and subsequently decoded in a step that utilizes an acquisition gradient applied 
while the signal from the spins is being captured in a recording step. The application of 
a proper acquisition gradient can successfully decode the information that was encoded 
during the course of the initial evolution time, in the form of an observable echo signal. 
The timing of this echo will depend on the strength of the Qi internal interaction that 
acted on the spins prior to the mixing step, hence allowing mapping of the spectrum 
along the indirect domain. Furthermore, this unwinding step can be immediately 
reversed, and then, repeated multiple (N2) times by alternating the sign of the 
acquisition gradient, thereby allowing monitoring of the D2 frequencies of the spins 
active during the second, directly-detected tg period. Signals obtained during such cyclic 
rephasing/dephasing train can be arranged into a bidimensional data set, which by 
Fourier analysis along t2 will lead to a desired 2D NMR spectrum correlating and Q2 
frequencies. The incorporation of multiple linearly-independent gradient geometries 
also enables an extension of this protocol to arbitrary N-dimensional NMR processing or 
experiments, where frequencies along N-1 indirect domains are spatially encoded using 
the new method and the final time-domain is monitored in the usual direct fashion. 
Fourier analysis of the spectral peaks along the indirect k/vi dimension provide as a 
fringe benefit the spatial positioning of the associated metabolite, thus yielding MRI 
information within the same experimental set. 

There are many applications and uses of the present invention. To render a more 
comprehensible understanding of the invention, the following description will be with 
particular reference to NMR. However, it should be clearly understood that the 
applicability of the invention is to any type of spectroscopy as enumerated in the 
foregoing. The ensuing discussion will be limited primarily to NMR for the sake of 
simplification and not by way of any limitation of the scope of the teachings herein. 

Recent developments in probehead and magnet technologies have increased 
NMR's sensitivity by almost an order of magnitude during the last decade, see R. F. 
Service, Science 279, 1127 (1998) and A. Constans, The Scientist 17, 45 (2003), thus 
limiting the acquisition times of many multidimensional NMR experiments to signal 
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digitization rather than to signal-to-noise (S/N) considerations. There are several areas 
where this bonus in signal enhancement is combinable with significant reductions in 
acquisition times brought about by the new method and apparatus of the invention, to 
speed up existing investigations or make new kinds of hitherto unfeasible investigations 
possible. The subordinate methods of the present invention include 

i. The novel method of the analysis by multidimensional NMR of rapidly-changing 

dynamic systems. The possibility of completing the collection of 
multidimensional NMR spectra within a 0.1 sec time scale will enable our 
invention to monitor in real time, a variety of chemical and physical processes 
and reactions that are hitherto outside the capabilities of NMR. These include 
the real-time monitoring of ongoing chemical reactions, and the folding of 
biological macromolecules. 

ii. The novel method for the application of multidimensional NMR to hyperpolarlzed 

spin states. A number of methods have been developed to enable the 
generation of very highly polarized spin states. These systems can impart on 
atoms and molecules NMR signals that are ca. five orders of magnitude more 
intense than a conventional NMR signal, but they are transient states that 
decay relatively rapidly and take long times to be generated. Only single-scan 
experiments are thus usually implemented on such hyperpolarlzed spin 
systems. The present invention enables the routine application of complex 
multidimensional NMR experiments to such hyperpolarized systems, enabling 
extensions of chemical studies. 

iii. The novel method for the characterization of analytes subject to flow through a 

NMR spectrometer, and thereby the coupling of multidimensional NMR with 
high-throughput chromatographic techniques. The combination of NMR with 
chromatographic techniques opens one of the most promising routes to the 
characterization of chemical and biochemical samples. The residence time of 
such flowing samples through the NMR reception coil, however, is very limited 
(«1 sec). Therefore only unidimensional NMR spectra have been so far 
collected in real time on this kind of samples. The method of the present 
invention will enable the acquisition of multidimensional NMR spectra on 
samples being chromatographed, thereby providing a new and much more 
powerful way to characterize plant extracts, natural products, amino acids, 
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peptides, nucleic acids and other types of chemicals being separated in a 
chromatographic column, 
iv. The novel method for rapid survey of large numbers of chemicals, like those 
made nowadays available by Combinatorial Chemistry. Combinatorial 
Chemistry is a novel approach to the synthesis of organic, inorganic and 
pharmacological molecules, whereby thousands of compounds are 
synthesized and tested in a variety of ways for chemical and/or biological 
activity. Combinatorial methods have provided much of the impetus for the 
ongoing revolution currently undergoing in Proteomics and Metabonomics. 
The enormous number of compounds that this approach requires be tested 
only allows high-throughput analytical techniques to participate in these tests 
and characterizations. The invention described herein allows the 
incorporation of ultrafast multidimensional NMR methods to this array of high- 
throughput techniques, thereby providing a new route to the discovery of new 
catalysts, new pharmaceuticals, pharmaceutically-active peptides and nucleic 
acids, etc. 

V. The novel method for the acceleration of quantum computing algorithms. NMR 
offers one of the most practical approaches to implement a quantum 
computer. The present invention enables the speeding up of such 
multidimensional NMR quantum computers by several orders of magnitude. 

vi. The novel method for the structural elucidations of large molecules, and in 

particular of systems involving proteins and nucleic acids. NMR-based 
structural elucidations on high molecular weight, complex systems will 
eventually demand the use of a large number of spectral dimensions (over 4) 
for achieving sufficient spectral resolution of the peaks. At the same time, 
however, large molecules in general and biological macromolecules under 
physiological conditions in particular, are usually incapable of withstanding the 
long acquisition times hitherto associated with such experiments. By 
speeding the times required to implement such experiments by several orders 
of magnitude the present invention bypasses such limitation, providing a new 
way to obtain the structure of macromolecules in their native states. 

vii. The novel method of the present invention for conducting in vivo spectroscopy 

and the following of fast metabolic processes using ultrafast multidimensional 
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NMR. Indeed, multidimensional NMR spectroscopy on animals and/or 
humans is currently hampered by the long times that subjects need to reside 
within the NMR magnet for the completion of the experiments, a demand 
which should be greatly eased by the practice of this novel invention. This will 
enable a new route to the clinical diagnosis of disease whether as a pure 
spectroscopic tool or in combination with methods for spatial localization. 

viii. The novel method for conducting MRI protocols, relying on magnetic field 

gradients to both speed up the acquisition of the spectral data as well as for 
locating the spatial position of spins in a non-invasive fashion. This opens new 
routes for the accelerated acquisition of clinical diagnostic and research MR 
images. The new ultrafast multidimensional MRI method of the present 
invention that results can be employed to monitor brain metabolism, pulsating 
regions (thorax, abdomen), etc. It can also aid for the real-time positioning of 
malignancies and hence as aid in surgical procedures. 

ix. The novel method to extract from multidimensional spectra collected within a 

single continuous scan, information about the spatial localization of the 
sample using a post-processing of the collected data. Multidimensional 
spatial localization can thus be combined with multidimensional spectral 
information within the same short time scale. This will enable a new approach 
to clinical diagnosis and basic physiological research. 

x. The novel method to shorten by orders of magnitude the duration of 

unidimensional NMR spectra collected for iow-y nuclei, by monitoring their 
shape indirectly via their more sensitive neighbors. For instance when 
dealing with heteronuclear ^^N NMR studies of organic molecules, the 
sensitivity of single-scan experiments based on using the invention for 
monitoring their 1D NMR spectra via the signals of neighboring nuclei 
could increase by ca. 5-10. This should greatly facilitate the research of 
medicinal chemists, utilizing heteronuclear NMR as a main tool for the 
preparation and characterization of pharmacologically-active and natural 
products. 

xi. The novel method to enable the acquisition of high resolution NMR spectra even 

when dealing with inhomogeneous or unstable magnetic fields. A 
considerable cost of current NMR and MRI instrumentation, stems from the 
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need of building intense magnets possessing homogeneities in the order of 1 
part-per-million (ppm) over the specimen volume. Only thus can high- 
resolution NMR spectra be resolved. Furthermore these magnetic fields need 
to remain unchanged ("locked") over the course of the whole spectral data 
acquisition, least their change ends up broadening the spectral response from 
the spins. By departing from the usual way of encoding and collecting NMR 
data the invention can alleviate both of these demands. High resolution 
spectra even when dealing with inhomogeneities that are orders-of-magnitude 
larger than those currently tolerable. This can have important consequences 
in the hardware used to collect uni- and multi-dimensional NMR spectra: 
much less sophisticated "tabletop" magnets could become routine for 
analytical use in high resolution NMR characterizations, and electromagnets 
whose stability and homogeneity are currently too poor (including the 
extremely high field Bitter magnets) could become suitable for the much 
higher demands of NMR and MRI 
xii. The novel method for extending single-scan multidimensional spectroscopy to 
techniques other than NMR, such as, mass spectrometry, electron 
paramagnetic resonance, as well as a variety of pulsed infrared, optical and 
UV spectroscopies. As in the case of NMR, these techniques are based on 
monitoring the response of the system as a function of an incremented time 
variable, thereby requiring the collection and processing of numerous 
individual scans. By adapting the ideas of a reversible inhomogeneous 
frequency broadening coupled to the selective excitation of sub-ensembles 
within the sample, the novel method of the present invention enables other 
kinds of multidimensional spectroscopies to be reduced to a single-scan 
acquisition. 

Other and further objects and advantages of the present invention will become 
more readily apparent from the following detailed description of preferred embodiments 
of the invention when taken in conjunction with the appended drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Figures 1A-B present graphically and pictorially steps of the invention as they 
relate to the underlying the acquisition of 2D NMR spectra within a single scan. 

Figures 2A-B present graphically and pictorially the steps of the invention as they 
relate to carrying out the selective excitation step in the initial stage of Figure 1, without 
encoding excitation artifacts, by means of a gradient-echo refocusing. 

Figures 3A-C present graphically and pictorially the steps of the invention as they 
relate to the acquisition while in the presence of an oscillating gradient; positioning of 
data points in the (k/vi,t2)-space; and further repositioning of the collected data points in 
preparation for a fast Fourier transform processing. 

Figure 4 is a simplified pictorial diagram describing the origin of peaks along the 
indirect dimensions of ultrafast NMR experiments. 

Figures 5A-G show graphically and pictorially the data acquisition and processing 
stages involved in the retrieval of ultrafast NMR experiments on a chemical sample. 

Figures 6A-D compare graphically and pictorially the similar results afforded by 

conventional vis-a-vis ultrafast 2D ^H-^H COSY and TOCSY NMR on a chemical 
sample. 

Figures 7A-D are similar to Figure 6, but for cases involving two different kinds of 
heteronuclear 2D correlations: (A) with direct detection; (B) with indirect detection. 

Figures 8A-B illustrate the inventive method's capability to the implementation of 
constant-time 2D NMR experiments, using simple pulse sequences as examples. 

Figures 9A-E illustrate the applicability of the inventive method to the acquisition 
of 2D NMR images (pure 2D MRI). 
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Figure 10 illustrates the possibility of the invention to retrieve both spectral and 
spatial multidimensional information, from the same single-scan sub-second acquisition 
(spatially-resolved 2D NMR/MRI) 

Figure 1 1 illustrates conceptually the invention as shown in Figure 1 regarding 
the application of the invention concerning the ultrafast acquisition of 3D NMR spectra. 

Figures 12A-D demonstrate the inventive method illustrated in Figure 11, with the 
single-scan acquisition (M40 ms) of a 3D NMR spectrum on glycerol/D20. 

Figures 13A-B shows the present invention as applied to the acquisition of 
arbitrary A/-dimensional NMR spectra (A), and illustrates that according to the invention 
a single-scan 4D NMR experiment can be completed within 94 ms (B). 

Figures 14A-C demonstrate the potential of the invention to increase the 
sensitivity of unidimensional NMR acquisitions on heteronuclei (i.e., nuclei other than 
protons) by relying on spatially-selective indirect-detection protocols. 

Figures 15A-C give a schematic description of how the invention could be used to 
retrieve high resolution NMR spectra even when dealing with inhomogeneous magnetic 
field, by manipulating the phases of the RF pulses used in the slice selection for 
implementing a compensation in the field strengths experienced by spins in inequivalent 
slices. 

Figure 16 is a block diagram of a preferred embodiment of the apparatus of the 
present invention with respect to NMR/MRI. 

Figure 17 is a flow chart illustrating the inventive method as it is applied towards 
the acquisition of multidimensional NMR spectra. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE INVENTION 

Referring now to the drawings, preferred embodiments of the present invention 
will be described in detail. 
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The paradigm for 2D NMR as laid down by the core contributions of Jeener and 
Ernst, see J. Jeener, in Ampere International Summer School II, Basko Polje, 
Yugoslavia (1971) and W. P. Aue, E. Bartholdi and R. R. Ernst, J. Chem. Phys. 64, 
2229 (1976), is sumnnarized by the well-known scheme 

Preparation/Excitation - Evolution (ti) - Mixing - Detection (t2) (1 ) 

However, according to the present invention, a 2D spectrum I(vx,v2)can be 
obtained if this approach is followed to collect a two-dimensional signal S as a function 
of ti and f2. and subsequently Fourier transformed. Yet the signals that are actually 
generated by the nuclear spin ensemble can only be directly detected as a function of 
f2, and hence the spins' coherent evolution as a function of the remaining time variable 
cannot be monitored directly. According to the invention, the remaining time variable is 
monitored indirectly in a unique manner. Independent acquisitions are thus carried out 
where the value of ti is systematically incremented; since the mixing process step 
occurring prior to t2 preserves the encoding imparted during the ti evolution, Fourier 
transformation (FT) of the digitized signals as a function of the ti parameter reveals the 
nature of the Qi frequencies that affected spins prior to the mixing. 

It is clear that this 2D NMR scheme, involving homogeneous and sequential 
evolution, mixing and acquisition periods for all spins throughout the sample, is not 
amenable to implementation within a single scan. The ultrafast 2D NMR methodology 
that is proposed by the present invention bypasses this limitation by partitioning the 
sample into a series of independent subensembles, each of them characterized by an 
individual ti evolution. Such partitioning can be done, for instance, by applying a time- 
incremented series of Ni selective excitation sequences throughout an 
inhomogeneously-broadened sample (Figure 1). A homogeneous mixing period (Figure 
1 A) followed by the observation of the independent signals generated by each member 
of this set, can then enable the simultaneous acquisition of a complete conventional- 
type 2D NMR data set (Figure 18). The overall acquisition scheme would then be 
summarized as 
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r Selective 1 Inhomogeneous Homogeneous Inhomogeneous 

[Preparation/Excitation — Evolution{t{) ~ Mixing ~ Acquisition{t2) 

Though evidently more complex than the original 2D NMR scheme in eq. (1), this 
partitioning of the sample offers the potential to collect the complete bidimensional data 
set within a single scan. This of course, depends on techniques for endowing the 
various A/-/ sub-ensembles with individual ti evolution times, and then for monitoring the 
signals that these generate independently and within a single transient. 

Figures 1A and 1B show graphically and pictorially the method of the present 
invention and more particularly, depict a simplified scheme or method respecting the 
application of spatially-heterogeneous excitation and detection schemes, toward the 
acquisition of multidimensional NMR data within a single scan. The excitation that 
triggers the initial spin evolution is assumed to affect spins in different positions within 
the sample at a series of incremented times (Figure 1A). This creates a spatial 
encoding of the initial f ^ evolution period, which is then monitored as a function of t2 via 
the spatially-resolved acquisition of the NMR signal. Following a suitable processing of 
the signals monitored as a function of the t2 and spatial variables, this strategy allows 
one to retrieve a complete 2D NMR spectrum within a single scan (Figure 1B). 

It is possible to propose several alternatives for achieving the spatially 
heterogeneous evolution and detection processes described in eq. (2) and Figure 1. 
Out of the possible alternatives, the invention concentrates for simplicity on those 
illustrated in Figures 2A and 2B for the excitation process, and in Figures 3A, 38 and 3C 
for the acquisition stage. Figure 2A shows that a train of frequency-shifted RF excitation 
pulses possessing different offsets 0/ can be applied while in the presence of a 
magnetic field gradient, in order to achieve the incremented evolution (ff) of spins 
throughout the sample assumed by the protocol in Figure 1A. The train of frequency- 
shifted pulses needs to be applied in combination with a synchronized reversal in the 
sign of the field gradient (gradient echoes) if an evolution that is solely dictated by the 
internal spin evolution frequencies is to be achieved. Assuming that gradients in this pair 
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are applied for equal time lengths Tp (e.g., the RF pulse duration), this will result in a 
gradient echo where any dephasing that may have affected hitherto excited 
magnetizations becomes compensated. 

Figures 3A, 3B and 3C illustrate the protocol employed in order to achieve the 
spatially-resolved data acquisition stage that is required by the method or scheme as 
given in Figure 1A. Dots symbolize the data points digitized during the course of 
positive and negative acquisition gradients. This oscillatory gradient module is then 
repeated A/2 times, with A/2 defining the number of effective points along tz. In Figure 3B 
data points monitored during the hybrid spectroscopic-imaging module defined in Figure 
1A need then to be separated according to their respective k and t2 variables, with 

k= y^G„{t)dt . The open and closed dots symbolize the coordinates of such points. 

In Figure 3C, as data stemming from the rearrangement illustrated in Figure 3B are not 
arrayed within a regular grid ready to be processed using the fast FT algorithm, points 
are sorted once again into two independent bidimensional data sets that are then ready 
to be individually processed. Both 2D data sets can then be co-added after their 
processing for the sake of improving the overall S/N of the method or experiment. 

According to the invention, when dealing with a conventional isotropic sample, 
the simplest route to achieve an inhomogeneous set of evolution times is by imposing 
an auxiliary gradient on top of the homogeneous external magnetic field Bq. This in turn 
enables a sequential excitation of the spins, based on the use of a train of frequency- 
shifted selective RF pulses (Figure 2A). Using a gradient-based step or scheme 
requires imparting on successively-excited spin packets an evolution phase that reflects 
solely the internal coupling frequencies, but not the frequencies defined by the artificial 
gradient. This step or goal is achieved by the invention by following the action of each 
selective pulse with reversal of the +Ge gradient employed to implement the spatially- 
heterogeneous excitation via an opposite gradient of amplitude -Ge (Figure 2B). 
Assuming that gradients in this pair are applied for equal time lengths Tp (e.g., the RF 
pulse duration), this will result in a gradient echo where any dephasing that may have 
affected hitherto excited magnetizations becomes compensated. Alternatively a 



13 



continuous, frequency-swept (chirped) excitation may be used to achieve a similar goal 
of a encoding along a spatial dimension the time incrementation of the spins. 

The next aim is to monitor the signals originating from these individually-excited 
slices, after they have been subject to the mixing process step. This spatially-resolved 
detection of the signal as a function of t2 amounts to a hybrid spectroscopic/imaging 
experiment, an acquisition that can be implemented within a single scan using a number 
of alternatives, see for example, P. T. Callaghan, "Principles of Nuclear Magnetic 
Resonance Microscopy" Oxford University Press, Oxford, 1991, and B. Blumich, "NMR 
Imaging of Materials" Oxford University Press, Oxford, 2000. In particular, during this 
stage of the methodology's development, the collection of NMR signals while subject to 
an alternating field gradient (see Figures 3A, 3B and 30) was adopted towards this end, 
see for example P. Mansfield, Magn, Reson. Med. 1, 370 (1984). This leads to signals 
becoming simultaneously encoded according to the spins' spatial positions, as well as, 
according to their internal evolution frequencies. Such encoding can be summarized by 

an evolution phase <i>= ' where k= 1^ y^G^(t)dt is a gradient-related wave 
number encoding position, and t2 reflects the extent of the free evolution and encodes 
the internal frequency. Evolution while in the presence of an oscillating ±Ga gradient can 
then be represented by a "zig-zagging" trajectory throughout the (k,t2)-space, which is to 
be digitized with a sufficiently short dwell time to characterize both the /c- and f2-behavior 
in a continuous fashion (Figures 3A, 3B). Dots in such Figures symbolize the data points 
digitized during the course of positive and negative acquisition gradients. This oscillatory 
gradient module is then repeated N2 times, with A/2 defining the number of effective 
points along tz- Processing such set into a spatially-resolved 2D NMR spectrum requires 
discriminating the signal contained in these points against the two extraction variables 
(k,t2)' These points do not appear distributed over the regular 2D grid that is needed for 
carrying out a fast FT along f2, a complication that can be resolved by implementing a 
separation procedure like the one illustrated in Figures 3B, 30 prior to the transform. 
Separating interleaved data sets in such fashion provides two S(k,t2) signals displaying a 
regular spacing between consecutive points, thus being ready for the fast FT. The 
particular example detailed in the Figure assumes an oscillatory gradient with a square 
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shape; the actual shape, however, is immaterial, and certain instances may arise where 
alternative shapes (such as sinusoidal modulations) may become more advantageous. 

Figure 4 is a simplified pictorial diagram describing the origin of peaks along the 
indirect dimensions of ultrafast NMR experiments. The heterogeneous nature of the ti 

evolution leads to an encoding of the internal precession frequency Oi along the z axis 

(second panel from left); this spiral of spin-packets is subsequently unwound by an 
acquisition gradient Ga possessing an identical z spatial dependence. The coherent 
addition of spin-packets thus leads to a sharp echo along the k coordinate whose 
position reveals the extent of Oi encoding prior to the mixing process — in essence, the 
spectrum along the indirect dimension. Such "peak" formation is only illustrated here for 
a portion of the first acquisition gradient echo; the phase encoding gained by this echo 
peak during the course of the A/2 gradient-reversals occurring as a function of f2. 

provides a conventional route to measure the /22 frequencies active during the 
acciuisition. 

An odd feature resulting from this 2D NMR data acquisition and processing 
modes or steps, is the fact that the k axis automatically becomes the vi domain that 

defines the spins' indirect evolution frequencies. Thus by contrast to conventional 2D 
NMR, which requires FT against the ti parameter to characterize the Qi frequencies in 

the spin system, no FT of the data along the /c-axis is needed altogether This non- 
conventional processing stems from the fact that in the present invention and in this kind 
of experiments, peak positions along the indirect domain are not defined via the Fourier 
analysis of a time-evolving coherence, but via the displacement observed in /c-space for 
the position of an echo created by interfering magnetizations from different slices. A 
graphical depiction of this aspect or feature is illustrated in Figure 4, which analyzes the 
fate of magnetizations throughout a portion of an ultrafast acquisition assuming a 
simplified "five-slices" sample. The initial excitation segment of the sequence affects 
successively the various slices and imposes on them a spatial encoding; yet this 
encoding will only reflect the initial Di internal evolution frequency, since the effects of 
the auxiliary Ge are being compensated by reliance on gradient echoes. The extent of 
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the ensuing evolution can therefore be described for a particular site by a "winding" of its 
magnetization through the sample, with a pitch dictated by exp[iCD.i(z-zo)]. This spatial 
winding is then preserved as either a phase- or an amplitude-modulation throughout the 
mixing period, with the result that the overall signal available for detection at the 
beginning of the acquisition will in general be null. At this point, however, a gradient is 

applied on the sample, whose spatial exp[iy^ lGa(f)dt'z] dependence is capable of 

0 

unwinding the initial spiral of magnetizations. Such unwinding will lead to a constructive 
interference among the spin-packets corresponding to different slices whenever 

k=r^lG^(t')dt'=-CQ^, thus decoding the initial evolution frequency that was active 

0 

during ti via the location of an echo along the /c-axis. This echo makes up the "peak" 
that will be observed along the indirect domain, whose further evolution as a function of 
the time t2 is then monitored by successive dephasing and rephasing processes driven 
by changes in the sign of Gq. Hence no additional /c-axis FT becomes necessary in 
order to decode the indirect evolution frequencies active during ti. 

Such arguments can be employed to derive a mathematical expression for the 
observable shapes observed for peak as derived according to the invention or in this 
kind of experiments. To do so, relaxation and Q2 shift effects are neglected in the 
beginning, and its own individual magnetization vector is associated to each of the Ni 
discrete slices excited by selective pulses during the course of the initial evolution 
period. According to the method of the present invention, these spin-packets will end up 
equally-spaced throughout the sample's length L at positions 

\zj = — — [(A^i~1)/2- m , and contributing with identical weights to the overall 

S(k) signal that is detected. When acted upon by the acquisition gradient Ga, the spin- 
packets within such discrete slices will be endowed with instantaneous precession 

frequencies vy = zjyaGa- The signal originated as a function of k= \G^{V)dV can then 

0 

be written as 
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Sik)= IV''^ . (3) 



The {Ffij=0,Nl-1 coefficients here are a complex set describing the state of the spin- 
packets at the /c = 0 instant; they are the x-y magnetizations associated to the various 
spatial slices, to which is ascribed a y-dependence which accounts for the Q1 frequency 
that may have been active during the initial ti evolution period. Such spin-packets can 

consequently be written as = — e^^^^^^^\ A being the overall magnetization 

expected from a particular site when considered over the whole sample. On the basis of 
Figures 1-3 one can further express the evolution times ti experienced by 

magnetizations within each of these slices as /i(zy)=C (zy-z_i), where the 

coefficient depends on the spacing between consecutive selective excitations. 
Assuming as done earlier an ideal dwell time Ati equal to twice the excitation pulse Tp 

as well as yeGe/AO > 0 leads to C= — ^ . Using all these definitions for zy, Fy, ti 

and C, it is possible to rewrite eq. (3) as 



«T (yv,+i)^ ^1 



wfiere s = %l{rp{Ni -\)+kL. The geometric series in eq. (4) can be evaluated using the 
well-l<nown relation TIjSqz' =(1-z^)/(1-2), leading to 



which can also be rewritten as 
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As the ratio in the middle of this expression is unity, eq. (6) yields 



^ (A^,+i)n, sinj^^eA^i /(^i - 1) 



sin -e/(A^i-l) 



(7) 



On considering H'\»1 we have A/f+t «A/t, Ni/(Ni-1) »1 and 1/(Ni-1) =1/Nr, using 
these approximations and replacing back for c's expression results in the relatively 
simple line shape expression 



S(k) = 



^JT,Nfl, sini 



^1 



sin 



(8) 



These equations also enable a derivation of basic Nyquist criteria for the new 
acquisition step or scheme. Denoting 7a as the duration during which ±Ga gradients 

are applied and At as the actual physical dwell time (Figure 3) leads to Nk = To/At 
complex points being collected along the /c-axis. Assuming then that the ±Ga cycle is 
repeated A/2 times, results in a 2TaN2 total digitization time for the step or experiment 
and in 2N2Nk complex points being collected overall. Points following the data 
separation procedure will be spaced along t2 by 2Ta intervals, leading to SW2 = (2TJ~^ 
spectral widths along the direct domain. The maximum range that under such conditions 
will be scanned along the /c-axis is kmax = raGaTa, which will in turn enable peaks to be 
accommodated over a SW^ =|^max/Cl spread along the indirect domain. The constant C 
is given by the ratio between the extent of ti evolution and the degree of spatial z 
encoding: C=^&t^/Az = At^ /[AO //jG^], At^ being the time between one excitation pulse 
and the next, Ye^e gradient's strength during excitation, and AO = |0/+f-0/| the 
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AO 



These expressions for SW^ and SW^ are of practical quantitative 



constant increment used to offset the selective excitation pulses. Under the idealized 
scheme in Figures 1-3 A/, w2J^, twice the selective pulses' duration; C is then 

2TpYeGe^^0 and the overall indirect spectral width becomes (in frequency units) 

use, provided that corrections are made for the gradients' non-idealities. These 
arguments also predict that, at least under ideal conditions, purely absorptive line 
shapes are built-in features of ultrafast NMR experiments. Indeed peaks along the 
indirect dimension form as a result of the constructive interference among spin-packets 
distributed along the rectangular p(z) profile of the sample container. Leaving potential 
distortions arising from the selective excitation and the spin relaxation aside, this will 
lead to peaks possessing sinc-type point-spread functions and thus no dispersive 

component. When convoluted with the normal Absorptive(v2)-^i'Dispersive(v2) line 

shape resulting along the direct domain after FT(f2), a 2D sinc(vi)-Absorptive(v2) kernel 
possessing convenient site-resolving characteristics will characterize peaks in this kind 
of spectra. 



Figures 5A to 5G illustrate and show graphically and pictorially the summary of 
events involved in the single-scan acquisition of phase-sensitive 2D NMR spectra, 

illustrated with data recorded on a solution of n-butylchloride dissolved in CDCI3 and 
utilizing a 2D NMR sequence devoid from actual mixing process. In Figure 5A time- 
domain data is collected using the spatially-selective excitation/detection procedures 
illustrated in Figures 2B and 3A. In Figure 5B, the magnified inset shows the signal 
(magnitude) arising from an individual Ta period, depicting in essence the compound's 

unidimensional vi spectrum. In Figure 5C, a 2D contour plot of the unidimensional data 
set illustrated in Figure 5B, following a rearrangement of its 2NkN2 points according to 
their k and t2 coordinates according to the procedure given in Figure SB. Interleaved 
data sets acquired with +Ga and -Ga gradients are still present at this point, thus 
resulting in a mirror-imaging of the signal along the /c-axis. Figures 5D and 5E show 
pairs of data sets resulting upon separating the interleaved +Ga/-Ga arrays in (C) into 
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two (k,t2) signals possessing NkN2 points each, as illustrated in the process leading 
from Figure 3B to Figure 3C. The signals shown in these sets have been subject to a 
phase correction and to a minor shearing that compensates for non-idealities in the 
acquisition gradient strengths (see below). The spectral structure observed already at 
this point along the k axis Is noticeable. Figures 5F and 5G show mirror-imaged 2D 
NMR spectra arising upon subjecting the data sets in Figures 5D and 5E to t2 Fourier 

transformation. The /c-axes of Figures 5D and 5E automatically become the axes of 
these spectra. 

As described above, Figures 5A to 5 G illustrates a basic 2D *H NMR acquisition 
on an n-butylchloride sample, using a 2D test pulse sequence (Figure 5A) that actually 
involves no mixing process or step as an example of the data processing required by the 
methodology of the present invention. Identical frequencies have here been active 
throughout the evolution and acquisition times, and only diagonal peaks result. As in all 
ultrafast experiments that will be later reported these data were collected as an 
unidimensional set, which when magnified appears to be composed by a series of 1 D 
indirect-domain NMR spectra (Figure 5B). This is a reflection of the method's built-in 
capability to FT the signal that was encoded during ti, via the Ga gradient defining the 
/c-axis. Processing continues with the assembling of the S(k,t2) 2D data set (Figure 5C) 
followed by the separation of the points corresponding to alternate +Ga/-Ga 
acquisitions, whose onwards processing is carried out independently (Figures 5D, 5E). 
A recursive artifact that may arise in ultrafast 2D NMR data sets consists of a "tilting" of 
the /c-space peaks, when followed as a function of time f2. This artifact grows as the 
intensity and/or duration of the gradients used in the experiment increases, and 
numerical simulations further detailed below suggest that it arises from a small 
weakening in the gradients' absolute value as a function of the acquisition time t2 
("gradient drooping"). Rather than compensating for this effect instrumentally, it can be 
corrected via a numerical shearing of the data. If purely absorptive line shapes are 
desired, processing continues after this stage by the phasing of the t2 - 0 slice, a 
copying of such phase correction through all values, and a Fourier transformation of 
the data as a function of to yield the final 2D spectrum (Figures 5F, 5G). Though not 
implemented for this data set, weighting and zero-filling can also be applied along the /c- 
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axis in a procedure that includes ancillary transformations of the data to and from their 
conjugate z-domain. As a result of all this processing, one obtains two sets of 2D (k, v2) 

data; essentially two NMR spectra, which are mirror-imaged to one another as a result 
of the opposing gradients employed in the scanning of their /c-axes. It is thus possible to 
Improve the overall S/N by reversing one of these sets and then combining them both 
into a single 2D NMR spectrum. As the and gradients defining these data sets 
are rarely of identical magnitudes, care should be exercised in this procedure and the 
coincidence of peaks in the two sets should be ensured prior to their recombination. 

The 2D NMR data in Figures 5A to 5G arise from a sequence deprived of a real 
mixing process, and hence its peaks are arranged solely along the main homonuclear 
diagonal. Figures 6A to 6D compare graphically and pictorially between single-scan 
(Figure 6A and 6B) and conventional (Figures 6C and 6D) 2D phased data sets 
acquired on an n-butylchloride/CDCl3 sample. Schematics of the COSY and TOCSY 
sequences utilized in these experiments are shown on the left- and right-hand side 
panels, respectively. The structure of the molecule, indicating its proton sites* shifts (in 
ppm), is also shown on top for ease of analysis. All data were acquired on a Bruker® 

DMX500 NMR spectrometer using a Nalorac® TXI probehead, and for each case the 
experiment's approximate duration is indicated. For the single-scan acquisitions Ni = 40 
initial Gaussian pulses were applied at 4 kHz offset increments while in the presence of 
YeGe = 150 kHz/cm, while the acquisition involved N2 = 256 gradient echoes with Ta = 
340 |js and 10 ps dwell times. All remaining pulses were applied non-selectively. 
Conventional 2D NMR experiments involved a TPPI phase cycle with 16 scans / ft point 
and 8 dummy scans. 

The comparison shown in Figure 6 illustrates the contrast between conventional 
and ultrafast 2D NMR spectra for the model compound when subject to two "real" 2D 
NMR sequences: COSY, correlation spectroscopy, which correlates cross-peaks among 
directly coupled neighbors, and TOCSY, total correlation spectroscopy, which 
establishes cross peaks among the full system of mutually-coupled spins, see for 
example R. R. Ernst, G. Bodenhausen and A. Wokaun, "Principles of Nuclear Magnetic 
Resonance in One and Two Dimensions" Clarendon, Oxford, 1987, and H. Kessler, M, 
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Gehrke and C. Griesinger, Angew, Chem, Int Ed, EngL 27, 490 (1988). As can be 
appreciated from Figures 6A to 6D, both normal and fast acquisition schemes convey an 
identical spectral information, though in radically different amounts of experimental time. 
The spatial encoding principles illustrated so far for monitoring 2D homonuclear 
connectivities, can also be employed for implementing heteronuclear correlations. 
Application of the invention in this respect or such experiments will rely on applying the 
Qf -driven spatial winding of magnetizations on a particular spin species S, preserving 
this encoding as an amplitude modulation during a mixing which implements a 
heteronuclear S^l transfer, and then decoding this information during the acquisition by 
gradients applied on the second / species. 

Figures 7A. 7B, 7C and 7D show examples where 2D correlations are 

established between different kind of nuclei. Figures 7A and 7B show a 2D "IH-I^c 
single-scan NMR spectrum obtained on a pyridine/CDCIs sample using the directly- 
detected heteronuclear correlation sequence sketched on top, which involves a spatial 
encoding of the magnetization and subsequent decoding on the ^^C channel. Data 
acquisition parameters included Ni = 33, Nk = 29, A/2 = 360, AO = 4 kHz, yeGe = 38 
kHz/cm, yaGa = 50 kHz/cm, 256 ^s Gaussian excitation rf pulses, 5 jus dwell times and 
homospoil (H.S.) gradients tuned to select the HxCz->HzCx coherence transfer. The 
total acquisition time was in this case «200 ms; other details are as in Figure 6. In 

Figures 7C and 7D, a 2D ''H-ISc pure-phase NMR spectra obtained on the indicated 
sample at natural abundance, using the indirectly-detected heteronuclear single- 
quantum correlation (HSQC) type sequence sketched on top. Data acquisition details 
included Ni = 20, Nk = 32, N2 = 128, zlO = 4 kHz, /eGe = 40 kHz/cm, /aGa = 97 
kHz/cm; other details are the same as mentioned in conjunction with Figures 6A to 6D. 
Two scans with the selective '^^C pulse (and the receiver) phase-cycled by 180' were 
collected for the sake of removing residual '^H-'l^c background signals. 

Figure 7 illustrates how two sequences that form the basis of several directly- and 
indirectly-detected heteronuclear correlation experiments, see for example R. R. Ernst, 
G. Bodenhausen and A. Wokaun, "Principles of Nuclear Magnetic Resonance in One 
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and Two Dimensions" Clarendon, Oxford, 1987; H. Kessler, M. Gehrke and C. 
Griesinger, Angew. Chem. Int Ed, Engl. 27, 490 (1988); J. Cavanagh. W. J. Fairbrother, 
A. G. Palmer III and N. J. Skelton, "Protein NMR Spectroscopy: Principles and Practice" 
Academic Press, San Diego, 1996; and M. Levitt, "Spin Dynamics" John Wiley & Sons, 
New York, 2001, can be modified for their single-scan execution according to the 

present invention. The first of these sequences is based on a evolution followed 
by an INEPT-based transfer period; the fact that ti in this sequence is triggered by a 
simple 90° excitation, enables the introduction of a spatial-encoding in ways that are 
entirely parallel to the ones implemented in the homonuclear correlations (Figures 7A 
and 7B). The only major distinction with the conventional 2D sequence then arises from 
the inability of conventional decoupling sequences to achieve their aim while under the 
action of intense oscillating gradients. Such goals were consequently implemented 

during both the and free evolution periods, by inserting 180" pulses within the 
complementary RF channel during the course of the +G/-G transition delays. The same 
decoupling principle was incorporated in the second of the heteronuclear correlation 
sequences assayed, HSQC (Figures 7C and 7D). By contrast to its directly-detected 
counterpart this sequence begins with a multi-pulse refocused INEPT block that is not 
susceptible to a simple spatial encoding. Still a spatial winding of magnetizations could 
be incorporated by shaping the 90" pulse that triggers the HzCz^HzCx conversion, and 
with it the ti evolution, following the first INEPT block. 

In addition to a stepwise ti incrementation, a common NMR approach to encode 
indirect evolution frequencies relies on systematically varying the location of a 180° 
pulse through the duration of a constant evolution period x, see for example Ernst et al 
and Kessler et al cited above. Such pulse will refocus the linear evolution terms 
(chemical shifts and heteronuclear couplings) yet leave homonuclear couplings 
unaffected, thus providing an appealing option to achieve homonuclear decoupling. 
The spatial encoding scheme introduced according to this invention can also be 
incorporated into this so-called constant-time modality, as illustrated by the example 

shown in Figures 8A and 88. Figures 8A and 88 show a 2D pure-phase single-scan 
NMR spectra collected on a CDCI3 solution of butyric anhydride, using the constant- 
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time modality sequences schematized on top of each contour. Data acquisition details 
included Ni = 26, Nk = 65. N2 = 160, zlO = 4 kHz, yeGe = 77 kHz/cm, yaGa = 46 
kHz/cm, dwell time = 8 ^s, 400 ^s Gaussian 180' pulses. Figure 8A shows a mixing-less 
pulse sequence, showing the three nonequivalent proton sites in the molecule aligned 
along the main diagonal. Figure 8B is idem but incorporating a 50 ms long isotropic 
mixing period, leading to cross peaks among all mutually coupled protons (right). 

This ultrafast scheme achieves a spatial encoding not by incrementing the initial 
excitation pulse but by using a single 90' excitation, which is subsequently followed by a 
train of spatially-selective 180' refocusing pulses. An exp[iCli(v'2Cz)] shift-driven 
winding of magnetizations is then imposed on the sample, from which the frequency 
spectrum can be decoded via the /c-space protocol just as described before. As 
illustrated by the butyric anhydride example shown in Figures 8A and 8B purely 
absorptive line shapes are again a possibility in this kind of acquisition mode. 

A particular class of 2D NMR experiments or applications are those underlying 
the acquisition of 2D magnetic resonance imaging (MRI). 2D MRI is, arguably, the most 
widely executed kind of NMR experiment, see for example M. A. Brown and R. C. 
Semelka, "MRI: Basic Principles and Applications" New York, 1999, yet from a 
spectroscopic standpoint it is actually a particular case of the much wider world of 2D 
NMR, see for example Ernst et al cited above. The main difference between 
conventional 2D NMR and 2D MRI applications or experiments is that whereas in the 
former the interactions to be correlated are purely internal, for instance couplings or 
shifts, the latter tends to neglect these for the sake of monitoring an artificial interaction 
given by the application of an external gradient in the magnetic field. Field gradients are 
also integral constituents of the spin evolution periods of the ultrafast scheme in Figures 
1A and IB, yet particular care is taken in this sequence to remove their effects via 
systematic changes in their sign and the ensuing generation of gradient echoes devoid 
of imaging information. If, however, these gradients oscillations were to be removed, 
the same acquisition scheme would become useful for the ultrafast acquisition of 2D 
MRI sequences. 
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Figures 9A to 9E show the adaptation and exemplification of how the new 
inventive method can be employed to collect 2D NMR images (pure 2D MRI) within a 
single scan. Figure 9A shows the modifications implemented on the basic spectroscopy 
scheme given in Figure 5A, whereby the position of the spins becomes encoded by 
removing the gradient echo refocusing from the excitation process. The net result is a 
train of Ni frequency selective pulses, applied while in the presence of an intermittent 
gradient. Figure 9B shows a pulse sequence resulting upon taking the Gq "off' period in 
Figure 9A to zero; the consequence is a continuous, highly efficient RF pulse whose 
offset varies linearly between Oi and Oni. Figures 9C and 9D show 2D MRI images 
obtained using the pulse sequences illustrated in Figures 9A and 9B, and reflecting as 
contour plots the water location profile used in the phantom that was experimentally 
tested (shown in the center for the sake of comparison), see Figure 9E. 

Figure 9A illustrates the basic single-scan 2D MRI experiment that would result 
from implementing these changes, with Qi becoming now associated with a position- 
dependent frequency 6/^ and encoding a position-dependent frequency G/2 (with 
and r2 reflecting any pair of x, y or z orthogonal directions). Furthermore, as no need for 
waiting for a refocusing -Ge delay is now needed after each initial RF excitation pulse, 
the possibility arises of employing a continuous, windowless train of variable-frequency 
RF pulses over the course of ti. Figure 9B presents the pulse sequence resulting when 
taking this concept to the limit of very short pulse widths; a single chirp pulse where 
frequency offsets are continuously swept between O1 and Oni then results, capable of 

carrying out the spins' excitation in a very efficient and robust manner. As illustrated 
with the test images presented in Figures 9C and 9D, either of these schemes is 
capable of affording 2D MRI images in a single, ultrafast scan. In this respect the 
abilities of this scheme are comparable to those of echo-planar imaging (EPI), a 
fundamental MRI tool that is also capable of affording 2D images within a single scan, 
see for example P. Mansfield, J. Phys. C: Solid State Phys. 10, 55 (1977), and M. K. 
Stehling. R. Turner and P. Mansfield, Science 254, 43 (1991). The principles of 
operation underlying these two methods are, however, entirely different, with the new 
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spatial-encoding approach of the present invention as applied to 2D MRI deriving from a 
general methodology capable of speeding up any 2D NMR sequence. 

Besides the application just described to MRI the invention provides a second 
route to obtain localization insight, this time via the spatial distribution information that 
resides in the shape of the spectroscopic echo peaks that are observed along the 
indirect k domains. To analyze this localization feature we focus on a particular 
(/2t, O2) peak, and examine the one-dimensional line shape that it will exhibit along the 
k/vi'3x\s as a function of the displacement Ak from the k = -C/37 condition. Using the 
bracketed summation in eq. (3) as starting point, leads then to 

The e '^^^^1 factor in this equation arises from the non-coincidence assumed between 
the z = 0 sample coordinate and the ft = 0 condition, and if so desired could be 

corrected away via an appropriate first-order phase correction. Aside from this artificial 
factor one recognizes in the S(Ak) function a Fourier conjugate of the spins signal 
amplitudes throughout the sample's z profile. These amplitudes — proportional in turn to 
the spins' z density weighted by the various spin relaxation processes that might have 
been active throughout the pulse sequence— can therefore be extracted via a Fourier 
analysis of the signal's line shape as a function of Ak, For instance when dealing with a 
homogeneous cylindrical sample and in the absence of spin relaxation, S(Ak) will display 
the Sinc-type dependence derived from eq. (8); Fourier analysis of this point-spread 
function against Ak will naturally reflect the square sample profile that originated it. In 
more general cases, a similar ancillary FJ{Ak) procedure on each (jOi, Q2) resonance 

would allow one to obtain the spatial distribution profile of the spins that lead to the 
formation of that particular peak. A most appealing aspect of such an approach is that 
spatially-localized information can thus be achieved without really demanding any 
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additions or modifications to the original single-scan 2D NMR spectroscopy experiment. 
All that is needed is extracting a suitable complex data array defining the full extent of a 
chosen peak along the k/vi axis, and subjecting it to Fourier analysis. Figure 10 

illustrates this mode of operation, with a set single-scan 2D NMR results recorded on 
the polar compound L-cysteine hydrochloride (8 mg) dissolved In an aqueous phase, 
and coexisting with an organic solution of ethyl acetate (5 ^iL) dissolved in CCI4. An 

isotropic recoupling sequence was incorporated during the course of the mixing leading 
to the single-scan acquisition of a 2D TOCSY NMR spectrum for the two solutes. 
Cross-peaks in this spectrum were sufficiently resolved to enable a spatial 
characterization of the analytes' distributions by means of their individual FJ{Ak) line 
shapes; the resulting profiles clearly show the spatial separation of the two phases for 
all the diagonal- and cross-peaks chosen for analysis. 

Technical aspects worth discussing are the spatial range and resolution that will 
characterize this kind of spatially-resolved 2D spectral protocols. These will actually 
follow from the spectral width (Sl/V) characteristics of single-scan 2D NMR experiments; 
for the basic setup summarized in the preceding paragraphs, their values along the 



indirect and direct domains have been shown given by sw. = ao^^^^^^ 



, SW2 = (2Ta)-i 



respectively. The spatial characteristics that follow from the Fourier analysis of peaks 
along the fc/vt-domain can be then worked out from coupling these SM/ considerations 

to the spectral characteristics of the analyzed metabolites, which will in turn dictate the 
maximum range of k values that can be used in the FT of any given peak. Indeed 
assuming that a peak to be Fourier analyzed along the indirect /c/vt-domain is resolved 

from other peaks along the same axis by a frequency span, will enable data points to 



be transformed over a maximum range of Ak^^ = 



. For the typical NMR 



AO 

parameters expected in a microscopy setting and at moderate magnetic field strengths 
kHz, AO «5 kHz, Ge =^10 G/cm, Ati ==1 ms) this will lead to a spatial 

resolution in the order of 0.2 cm. The associated field of view will in turn be given by the 
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number of points sampled within each Tq interval; typical dwell time values lead then to 
spatial fields in the 2-4 cm range. The overall order of magnitude of these two spatial 
parameters can be expected to increase or decrease by a decade upon switching the 
experiment to whole-body or microimaging conditions respectively; in either case, this 
simple analysis reveals that there is an appropriate window of opportunity where both 
spectral and spatial information can be extracted from single-scan experiments. 

The present description has focused so far on the use of spatial encoding 
methods and protocols to accelerate the acquisition of various 2D NMR experiments. 
The invention employs the same principles for acquiring higher-dimensional NMR 
spectra within a single scan, a particularly important goal given the exponential increase 
that the duration of A/-dimensional NMR experiments exhibit with respect to A/. A 
relevant point to notice in order to proceed with the extension of the new ultrafast 
acquisition invention to higher dimensionalities, are the dissimilar roles that the 
gradient's strength and the gradient's geometry play in the collection of the data. As 
evident from the arguments above, the strengths Ge, Ga will define important range and 
resolution characteristics of the single-scan 2D NMR spectrum. The actual geometrical 
distribution of the gradient by contrast, is mostly responsible for relatively minor line 
shape characteristics related to the sharpness of the echo formation. Indeed most of 
the arguments described in the preceding discussion would remain equally valid if the 
gradient used to encode and decode the frequencies were assumed applied along 
the X (or y) rather than along the z direction. In either case, a spiral of spin-packets 
would be created by the selective excitation protocol illustrated in Figure 1, even if for a 
cylindrical profile the use of radial (xory) gradients would end up giving different relative 
weights A(r) to the various spin-packets into which the sample is partitioned. In fact the 
only demand that the performance of this kind of experiments places on a gradient's 
geometry, is that it should generate a spiral of spin-packets that is convoluted enough to 
result in a null magnetization when considered over the whole sample. In principle there 
are infinite gradient geometries capable of fulfilling such conditions, including among 
others those represented by the boundless {Yl,mO,^} series of real spherical harmonics 
of the kind present in modern high resolution NMR shim systems, see for example M. J. 
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E. Golay, Rev. Sci. Instr. 29, 313 (1958). Thus if assuming a generalized field gradient 
possessing Yi^m spatial dependence and G = dBo/d(Yi^m) strength, all arguments 
leading to eq. (3) can be repeated to derive a macroscopic magnetization that at the 
conclusion of the ti evolution period will be summarized as 



This equation represents a winding of spin-packets, arranged this time along the 
gradient's Yj^m geometry and possessing non-uniform A(x,y,z) weights. Such winding 
will once again lead to an overall zero magnetization when considering an chemical 
shift evolution and integration over the whole sample. The same Yi^m gradient used to 
achieve the spatial encoding in eq. (9), however, holds the key for unwinding this spiral 
of spin-packets, and consequently, for the generation of an observable echo during the 
course of the acquisition. In order to achieve this unwinding, spins should be allowed to 

accumulate a phase exp(ihYi^m), with /r= J Ga(t)dt = -CQ^ during the course of each 
acquisition dwell time At2^ Just as in the z example above, this will again allow one to 
map the evolution frequencies that had been active prior to the mixing process in the 
form of a constructive interference among individual spin-packets. Also as in the z-case 
above, this spatial helix of spin-packets can be wound and unwound numerous times by 
periodically reversing the Yj^m acquisition gradient, in a process that will encode the spin 
frequencies as a function of t2 and thereby enable the collection of 2D NMR spectra 
within a single scan. The actual shapes of the resulting /c-echoes — shapes which will in 
turn define the kind of peaks observed in ultrafast experiments along the indirect 
domain — will depend on the sample's and gradient's spatial dependencies. The echo 
positions on the other hand will be independent of these details, and solely reflect a 
site's given shift. As can be appreciated from this Figure numerical simulations 
reveal a behavior that overall is uniform, with minor line shape differences that can be 
rationalized in terms of the characteristics assumed for the gradients and the sample. 




(9) 
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It follows from these considerations that only if identical gradient geometries are 
employed to impose the initial spatial encoding and the subsequent decoding of the 
internal frequencies, will a successful unraveling of the indirect evolution frequencies 
result. By contrast, the application of an intermediate gradient that is orthogonal to the 
one employed in the initial excitation, will neither add nor detract from the unwinding 
demands placed on the acquisition gradient. This in turn opens up a route for extending 
the protocol introduced in the preceding paragraphs to include the encoding of an 
additional indirect dimension, if modified as illustrated in Figure 11. Figure 11 shows a 
basic scheme proposed for the single-scan collection of 3D NMR spectra. A train of 
frequency-shifted pulses is first applied to achieve the spatial encoding of spins 
throughout different axial positions in the sample (ft); this is followed by a second rf- 
driven encoding of the spin evolution during t2 along a linearly independent radial 
direction. Data are finally collected while in the presence of oscillating acquisition 
gradients, which decode the initial Qt, CI2 frequencies along the (kz,kx) axes while 
monitoring the spin evolution along a third (^3) time axis. The actual direction (or even 
linearity) of the gradients used in the encoding is not fundamental. 

Such scheme incorporates two separate gradients arranged along linearly 
independent geometries, which implement two consecutive spatial encodings of the spin 
evolution. Each one of these encoding processes, assumed here for simplicity to lie 
along the x and z directions, proceeds independently and along an outline similar to the 
one described previously for the single-axis 2D NMR method or experiment. The first of 
these gradients will thus induce an Q^f^-dependent winding of the spin-packets along the 
z direction, while for each one of these z slices the second gradient will generate an 
Q2f2-dependent encoding along the x axis. Because of the ensuing "double-winding" of 
spin-packets the overall bulk magnetization is again reduced to zero, and an acquisition 
process implemented on the resulting sample will be associated with a null initial signal. 
Moreover, by contrast to the 2D encoding case, no single-axis gradient would be able by 
itself of regenerating an observable signal from such an initial state. Only the 
simultaneous application of Gz, Gx acquisition gradients can succeed in aligning the 
spin-packets, wound as they now are along two orthogonal encoding directions. In order 
to find the conditions leading to such an alignment the inventive method implements 
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again a step or protocol where both x and z gradients are rapidly and independently 
oscillated. The digitized signal can then be regarded as a function of three independent 
variables: ^ oc J G^iOdt ; ^ oc J G^iOdt ; and a time t3 associated with the final free 
evolution frequency Q3. When the ^-wave numbers happen to fulfill the kz = -CzQ^ and 
kx = -Cx02 conditions, spin-packets throughout the sample will interfere constructively 
and result in an observable echo. FT of this sharp 2D echo signal along the remaining 
t3 dimension provides then a means or way for characterizing all (Q^ CI 2, Q3) evolution 
frequencies experienced by the spins throughout the course of the method or 
experiment, from data collected within a single scan. Since the mechanisms by which 
peaks are originated in this kind of 3D NMR experiments are identical to the ones that 
were earlier described for ultrafast 2D NMR acquisitions, most of the arguments used to 
derive Nyquist and line shape criteria along the indirect dimension of the latter can be 
directly extrapolated to the indirect dimensions of the former. 

Figures 12A to 12D demonstrate the inventive method illustrated and described 
with reference to Figure 11. Figure 12A shows a TOCSY-HSQC pulse sequence 
assayed to corroborate the possibility of collecting a 3D NMR spectrum within a single 
scan. Figure 12B depicts an iso-surface representation of the full 3D NMR spectrum 
acquired on a glycerol/D20 sample, using the pulse sequence indicated in Figure 12A. 
The total time required to carry out this experiment was ca. 141 ms. Figures 12C and 
12D illustrate 2D contours resulting from projecting the 3D NMR data against the 
remaining spectral axis, and illustrate the expected line shapes for the compound. 
Excitation parameters included A/^ = 11, A/2 = 17, all selective pulses 256-^s-long 
Gaussians applied at 4 kHz offset increments, a 160 kHz/cm z gradient, and a 100 
kHz/cm X gradient. Initial and final heteronuclear transfer periods were set at 1.75 and 
1.40 ms respectively. Acquisition parameters included A/3 = 6, N4 = 24, Tq = 60 ps, 1.5 
|js physical dwell times, maximum x and z gradient strengths of 88 kHz/cm and 9 
kHz/cm respectively. Such conditions yielded 88x6x24 points in the mixed (v2/kz, vt/Zcx, 
ts) domain, which were processed as into a 128x16x64-point magnitude spectrum. 
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Figures 12A to 12D gives experimental proof on the feasibility of implementing 

single-scan 3D NMR, utilizing a 1H-13c-1h TOCSY-HSQC experiment on a 
glycerol/D20 sample, see as example J. Cavanagh, W. J. Fairbrother, A. G. Palmer III 
and N. J. Skelton, "Protein NMR Spectroscopy: Principles and Practice" Academic 
Press, San Diego, 1996. This heteronuclear correlation relies on the coherent nature of 
heteronuclear S->l transfer, which will preserve the spatial encoding imparted by a train 
of RF pulses even while the nature of the nuclei that carry such encoding changes, in 
order to parallelize the manner by which the heteronuclear evolution is encoded. Figure 
12A illustrates a pulse sequence based on such premise; Figure 12B shows the 
experimental results obtained with such sequence. 

Having discussed how the principles of single-scan 2D NMR can be extended to 
a third dimension and having illustrated the method for accomplishing this objective, the 
method of the invention can be adapted and can be extended to an arbitrary number of 
dimensions. The principles underlying the invention, as explained in detail above, can be 
generalized to an arbitrary number of dimensions. Such generalization is made possible 
by the countless gradient geometries that can, in principle, be utilized for encoding the 
spin evolution, represented for instance by the spherical harmonic set {Yj^m} discussed 
earlier. When applied in combination with a train of frequency selective RF pulses, any 
of these eigen functions can be employed to wind an independent spiral of spin-packets, 
of the type summarized by the Fourier relation in eq. (9). The task of unraveling the spin 
evolution frequencies that act along various indirect dimensions becomes then 
equivalent to scanning simultaneously the multiple ki^m associated to these various 
gradients. Such a procedure can in principle always be carried out within a single 
continuous acquisition, provided that flexible enough gradients are available. 

Figures 13A and 13B show the present invention applied to the acquisition of 
arbitrary A/-dimensional NMR spectra. In Figure 13A the method step is shown, of the 
present invention, for the acquisition of arbitrary A/-dimensional NMR spectra, and 
involving A/-7 independent spatial encoding events prior to the last mixing event, and 
their simultaneous decoding during the course of the signal acquisition. Figures 13A 
and 13B Illustrate the method of the present invention with respect to the acquisition of 
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a 4D NMR technique or experiment within a single scan, that can be completed within 
94 ms. As no actual mixing sequence is here involved between any of the tj f, +t 
intervals, the two peaks of this compound show as expected along the main diagonal of 
the 4D spectral space (grayscale isosurfaces). Excitation parameters included 19 
pulses applied along x and y axes (500 ps long), and 25 pulses along the z direction 
(250 ps long). The acquisition yielded 6x6x12x6 points in the mixed (vi/ky, v2/kxs 
t4) domain, which were processed into the displayed 6x1 6x32x1 6point magnitude 
spectrum. 

The portrayal in Figures 13A and 13B summarizes the kind of method or scheme 
that would then be involved in these arbitrary A/-dimensional experiments; the graph in 
Figure 13B illustrates experimental 4D NMR results collected on the practice of the 
disclosed method on the basis of such principles, with the triple-axes spatial encoding 
required by such sequence implemented with the aid of x, y, z spectroscopy gradients 
applied in combination with a constant-time protocol. Collecting similar spectra using 
conventional means might have taken several hours, perhaps even days, of continuous 
spectrometer use. 

The present invention also provides a new method to significantly shorten the 
times needed to collect unidimensional NMR spectra recorded from so-called 
heteronuclear spins X — basically any NMR targets of investigation other than protons. 
ID NMR studies on ^^N and ^^C nuclei in particular, have a great potential for chemical 
structure determination and play irreplaceable roles in both modern industry and 
academia. Because of the greatly reduced magetogyric ratios y characterizing these 
nuclei in comparison to ^H, however, the signal arising from ^^N and ^^C after a standard 
NMR experiment is usually less sensitive than its counterpart by a factor of {ynfyxf' a 
factor Yh/yx reflecting the smaller magnetic fields that the heteronucleus X can generate, 
and another yn/yx factor consequence of the higher Larmor frequency in which for a 
given field strength the NMR acquisition takes place. During the development of 2D 
NMR as a tool for biomolecular investigations, it was realized that this penalty in signal- 
to-noise could be compensated if heteronuclear experiments were carried out in 
"inverse detection" mode; i.e., if the heteronuclear evolution were encoded along the 
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indirect domain of a 2D experiment, and these information then passed over to a 
neighboring on which the actual measurement would take place. This mode of 
operation, however, finds limited application toward the acquisition of 1D ''^N and ^^C 
NMR data due to the fact that the organic or pharmaceutical chemist is usually not 
interested in paying the orders-of-magnitude longer acquisition times that are inherent 
to 2D over 1D forms of spectroscopy. The invention hereby presented, however, can 
easily overcome this time limitation. The starting point for doing so is again the pulse 
sequence illustrated in Fig. 1, but modified so as to "suspend" the information encoded 
along the direct domain. This essentially requires collapsing the direct-domain ''H NMR 
spectral information using a train of composite refocusing n pulses. This collapsing of 
the 2D NMR information into a single heteronuclear dimension should enable gaining 
considerable factor in sensitivity when compared with an analogous 1 D heteronuclear 
acquisition. When considering that for a given signal-to-noise ratio NMR acquisition 
times decrease as the square of an experiment's sensitivity, the potential savings that 
might result from implementing these sensitivity-enhanced 1D modes of acquisition 
become evident. 

Figure 14 demonstrates the feasibility of achieving significant signal-to-noise 
improvements in this manner, using 1D NMR of the ""^N nuclei as example. The 
sequence assayed for this test (Figure 14A) consisted of an initial INEPT block whereby 
''H z magnetization (Hz) was coherently transformed into two-spin HzNz spin order; a 
frequency selective train of ^^N radiofrequency pulses which in combination with the 
invention hereby proposed encodes the spin evolution of nitrogen spins along the 
sample's spatial orientation; a second INEPT-like period where this nitrogen encoding is 
transferred back to protons via a two-spin order state; and a final acquisition period 
combining a mono-polar decoding gradient with a train of tt pulses that effectively 
freezes the information and allows for the acquisition of a large number of identical 
^^N spectra within a single scan. All these spectra are then co-added for the sake of 
additional improvements in signal-to-noise. The indirectly-detected NMR spectrum of a 
25 mM ^^N urea sample dissolved in de-DMSO was collected in this manner, and 
compared with the standard direct acquisition usually used in 1D NMR. Whereas in a 
single scan an experiment based on direct monitoring of the ''^N spins the NMR signal is 
barely discernible from the background noise (Figure 14B), an indirectly-detected 1D 
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single-scan experiment leads to the presence of an unambiguous signal (Figure 14C 
and much higher signal-to-noise ratios. 

The present invention also provides a method for acquiring high resolution NMR 
spectra showing narrow resonance line widths, even in the presence of large static 
inhomogeneities or dynamic fluctuations in the main magnetic field. Achieving highly 
uniform, constant magnetic fields, with inhomogeneities of up to a few parts-pre-billion 
(ppb) throughout the sample volume and drift rates of the same order of magnitude over 
the duration of an experiment, is one of the most important and challenging demands 
posed to developers of NMR and MR! instrumentation. Magnetic field inhomogeneities 
or instabilities larger than this size preclude the retrieval of chemical shift spectral 
information, and complicate the achievement of high spatial imaging resolution. A large 
fraction of the cost of magnetic resonance Instrumentation is thus commonly invested in 
achieving these extremely uniform magnetic fields, employing advanced wiring materials 
and highly complex designs that include the addition of a considerable number of 
sophisticated "shimming" magnetic field gradients onto the main magnetic field. In every 
instance this shimming procedure is costly and time consuming. Moreover important 
occasions arise when shimming to ppb homogeneity is simply unattainable, for instance 
when dealing with MR! investigations on tissues with a high degree of anisotropic 
magnetic susceptibility or with "inside-out" setups of the kind used when employing NMR 
for water logging and oil prospecting (G. Coates, NMR logging: Principles and 
Applications, 1999). By opening an alternative to the usual acquisition of NMR spectra 
by directly detecting the spins' time-domain signal, the new invention hereby discussed 
provides a route for compensating for these artificial inhomogeneities. To explain how 
this compensation would operate we treat first the case of static inhomogeneities, which 
endowes spins at different positions in the sample with artificially different Larmor 
precession frequencies. The starting point for the compensation scheme lies then in 
encoding the evolution frequencies of subdivided spin-packets according to their spatial 
position using the principles illustrated in Figs. 1-4. We assume now that in addition to 
the natural spin evolution frequencies of interest there will be an encoding stemming 
from the undesired field inhomogeneities, which of course will not be at all as significant 
within the individual slices as it is on the much larger sample as a whole. Still, 
precession frequencies will adopt different values when considered from slice to slice 
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due to the latter's inequivalent positions and —just as they do in conventionally-detected 
time-domain NMR — affect the sharply-defined echo peaks that would othenA/ise be 
recorded (Figures 15A, 15B). On the other hand these inhomogeneities are static and 
thus their value can be precisely measured, for instance by mapping the spatial 
distribution of the main static magnetic field upon first installing it. Once the 
dependence of the magnetic field inhomogeneity AB with position zy is accurately known, 
a compensation for its effects can be easily introduced into the scheme of Figure 1 by 
manipulating the phases of the RF pulses involved in the spatial selection as illustrated 
in Figure 150. Phase corrections <fy ^-AB(Zj)[Czj+TQ] are here introduce, which 
compensate for the effects of static distortions over the sample as a whole but leave the 
effects of the internal evolution frequencies unchanged. These can thus be measured 
as detailed in Figures 3 and 4 of the present application. The main factor limiting how 
large a field inhomogeneity can be accounted for in this manner becomes the relative 
size of the field gradient G e used to impart the spatially-selective excitation, as this 
defines the "fine-ness" of the spatial selection and in turn the inhomogeneity within the 
slice. Preliminary calculations indicate that for conventional 1D inhomogeneity profiles 
field dispersions in the order of 100s ppm can be accounted for in this manner, and 
provide high-resolution NMR spectra in the presence of magnetic field dispersions that 
are orders of magnitude larger than those which can be nowadays be tolerated. 
Furthermore, the extension of a unidimensional scheme like the one illustrated in Figure 
150 to a three-dimensional case based on the concurrent utilization of x-y-z orthogonal 
field gradient sets, is straightfonA^ard. 

In addition of enabling for the compensation of static, pre-mapped 
inhomogeneities, the scheme described in the preceding paragraph could be used to 
compensate for dynamic magnetic field instabilities, and thereby enable the acquisition 
of 2D NMR spectra in such rapidly-drifting magnetic systems. This is a field of emerging 
importance, as the highest fields currently available for performing magnetic resonance 
measurements are generated by hybrid superconducting/Bitter designs relying on high- 
current direct-current (DO) power supplies. Due to the impossibility to regulate such 
intense current supplies down to parts-per-billion perfection the magnetic fields that they 
generate drift —usually by several ppm/sec. Such drift rates are incompatible with the 
acquisition of high-resolution NMR spectra in general, and with 2D NMR spectroscopy in 
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particular. The invention hereby presented will help alleviate these problems on the 
basis of two accounts: i) by shortening the overall times needed for completing a 
multidimensional NMR experiment by several orders of magnitude, instabilities occurring 
over a > 1 sec time scale cease to be a concern, and ii) by opening up the possibility to 
compensate for field instabilities by manipulating the RF excitation phases in the manner 
put forward in Figure 15C, short term instabilities (<1 sec) can also be taken care of. 
The first of these features follows clearly from the detailed description give above. To 
deal with the second of these features, the noise-like instabilities arising from DC power 
supplies, we propose carrying out the phase corrections detailed in the previous 
paragraph but in a real-time manner. This is conceptually analogous to correcting a pre- 
mapped field fluctuation, with the exception that the phases (|)(Z]) involved in the field 
compensation would have to be computed at the time of an RF pulse application. This 
in turn would need a circuit for monitoring the instantaneous changes in magnetic field 
occurring since the beginning of the NMR acquisition, something which can be 
implemented utilizing an ancillary electronic circuitry akin to the one normally 
incorporated in the deuterium lock circuits of high-resolution NMR instruments (see for 
instance E. Fukushima & W. Roeder, "Experimental Pulsed NMR: A Nuts and Bolts 
Approach". Addison, 1981). Multidimensional NMR studies in unregulated hybrid or 
electromagnets, could thus become a reality. 

The present invention also contemplates a method for the analysis or monitoring 
by multidimensional NMR of rapidly changing dynamic systems. As it is possible by the 
present invention to complete the collection of multidimensional NMR spectra within a 
0.1 sec time scale, the method of the present invention can enable monitoring in real 
time, a variety of chemical and physical processes and reactions that are hitherto 
outside the capabilities of NMR. These include the real-time monitoring of ongoing 
chemical reactions, and the folding of biological macromolecules, as described in C. M. 
Dobson and P. J. More, Nat. Struct. Biol. 5, 504 (1998). According to the present 
invention this is accomplished by a method for the real-time monitoring of a chemical or 
physical process comprising the steps of: (1) conducting a chemical or physical 
processin real-time; (2) monitoring the on-going chemical or physical process in real 
time by repeatedly acquiring multidimensional nuclear magnetic resonance spectra 
within a single scan at preselected short time intervals; each single scan being carried 
out comprising the steps of the methods described herein previously in detail. For 
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example, the single scan can be carried out by: (3) partitioning a predetermined sample 
of the on-going chemical or physical process into a set of independent subensembles; 
(4) applying the single scan to the sample by exciting the set of subensembles by a 
time-incremented series of selective excitation sequences; (5) generating a signal from 
each subensemble; (6) homogenously mixing the generated signals; and (7) 
simultaneously acquiring a complete multidimensional nuclear magnetic resonance data 
set for each preselected short time interval. 

The application of multidimensional NMR to hyperpolarized spin states is 
described by B. M. Goodson, J. Magn. Reson. 155, 157 (2002) and P. J. Carson, C. R. 
Bowers, D. P. Weitekamp J. Am. Chem. Soc. 123, 11821 (2001). A number of methods 
have been developed to enable the generation of very highly polarized spin states. 
These systems can impart on atoms and molecules NMR signals that are ca. five orders 
of magnitude more intense than a conventional NMR signal, but they are transient 
states that decay relatively rapidly and take long times to be generated. Only single- 
scan experiments are thus usually implemented on such hyperpolarized spin systems. 
The present invention can enable the routine application of complex multidimensional 
NMR experiments to such hyperpolarized systems, enabling extensions of chemical 
studies as described by B. S. Duckett, S. A. Colebrooke, Encyclopedia of Nuclear 
Magnetic Resonance, 9, 598 (2002), and of clinical studies as described by M. S. 
Albert, C. D. Gates, B. Driehuys, W. Happer, B. Saam, C. S. Springer Jr., A. Wishnia, 
Nature, 370, 199 (1994). This is carried out by the present invention by the method of 
for the real-time monitoring of a hyperpolarized spin system comprising the steps of: (1 ) 
monitoring an on-going hyperpolarized spin system in real time by repeatedly acquiring 
multidimensional nuclear magnetic resonance spectra within a single scan at 
preselected short time intervals; each single scan being carried out comprising the steps 
of the methods described herein previously in detail. For example, the single scan can 
be carried out by: (2) partitioning a predetermined sample of the on-going 
hyperpolarized spin system into a set of independent subensembles; (4) applying the 
single scan to the sample by exciting the set of subensembles by a time-incremented 
series of selective excitation sequences; (5) generating a signal from each 
subensemble; (6) homogenously mixing the generated signals; and (7) simultaneously 
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acquiring a complete multidimensional nuclear magnetic resonance data set for each 
preselected short time interval. 

The present invention has particular application to the characterization of 
analytes subject to flow through a NMR spectrometer, and thereby the coupling of 
multidimensional NMR with high-throughput chromatographic techniques. As described 
in H. H. Liu, C. Felten, Q. F. Xue, B. L. Zhang, P. Jedrzejewski, B. Karger, M. E. Lacey, 
R. Subramanian, D. L. Olson, A. G. Webb and J. V. Sweedler, Chem. Rev. 99, 3133 
(1999) and K. Albert, On-Line Liquid Chromatograptiy-NMR and Related Techniques. 
(John Wiley & Sons Ltd., Chichester, 2002), the combination of NMR with 
chromatographic techniques would be one of the most promising routes to the 
characterization of chemical and biochemical samples, if it were possible to do so. The 
residence time of such flowing samples through the NMR reception coil, however, is 
very limited («1 sec). Therefore only unidimensional NMR spectra have been so far 
collected in real time on these kinds of samples. By the method of the present invention 
acquisition of multidimensional NMR spectra on samples being chromatographed can 
be accomplished, thereby providing a new and much more powerful way to characterize 
plant extracts, natural products, amino acids, peptides, nucleic acids and other types of 
chemicals being separated in a chromatographic column. This is accomplished by the 
present invention by the method for the real-time monitoring of one or more samples of 
a material undergoing a chromatographic technique comprising the steps of: (1) 
monitoring an on-going chromatographic technique in real time by repeatedly acquiring 
multidimensional nuclear magnetic resonance spectra within a single scan at 
preselected short time intervals; each single scan being carried out comprising the steps 
of the methods described herein previously in detail. For example, the single scan can 
be carried out by: (2) partitioning a predetermined sample of the on-going 
chromatographic technique into a set of independent subensembles; (3) applying the 
single scan to the sample by exciting the set of subensembles by a time-incremented 
series of selective excitation sequences; (4) generating a signal from each 
subensemble; (5) homogenously mixing the generated signals; and (6) simultaneously 
acquiring a complete multidimensional nuclear magnetic resonance data set for each 
preselected short time interval. 
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The rapid survey of large numbers of chemicals like those made nowadays 
available by Combinatorial Chemistry. As explained in K. C. Nicolau, R. Hanko, W. 
Hartwig, Handbook of Combinatorial Chemistry: Drugs, Catalysts, Materials. (Wiley-VCH 
Verlag, Weinheim, 2002), Combinatorial Chemistry is a novel approach to the synthesis 
of organic, inorganic and pharmacological molecules, whereby thousands of compounds 
are synthesized and tested in a variety of ways for chemical and/or biological activity. 
Combinatorial methods have provided much of the impetus for the ongoing revolution 
currently undergoing in Proteomics and Metabonomics. The enormous number of 
compounds that this approach requires be tested only allows high-throughput analytical 
techniques to participate in these tests and characterizations. The invention described 
in this patent will now allow the incorporation of ultrafast multidimensional NMR methods 
to this array of high-throughput techniques, thereby providing a new route to the 
discovery of new catalysts, new pharmaceuticals, pharmaceutically-active peptides and 
nucleic acids, etc. This is accomplished by the present invention by the method for the 
rapid and real-time monitoring of a combinatory chemistry involving a plurality of 
samples of a number of different materials comprising the steps of: (1) rapidly 
monitoring an on-going combinatorial chemistry technique involving a plurality of 
samples in real time by repeatedly successively acquiring multidimensional nuclear 
magnetic resonance spectra within a single scan of successive samples; each single 
scan being carried out with respect to a sample and comprising the steps of the 
methods described herein previously in detail. For example, the single scan can be 
carried out by: (2) partitioning each successive sample of the on-going combinatorial 
chemistry technique into a set of independent subensembles; (3) applying the single 
scan to the sample by exciting the set of subensembles by a time-incremented series of 
selective excitation sequences; (4) generating a signal from each subensemble; (5) 
homogenously mixing the generated signals; and (6) simultaneously acquiring a 
complete multidimensional nuclear magnetic resonance data set for each successive 
sample. 

The acceleration of quantum computing algorithms is described by D. G. Cory, A. 
F. Fahmy, T. F. Havel. Proc. Natl. Acad. Sci. USA, 94, 1634 (1997). NMR offers one of 
the most practical approaches to implement a quantum computer. A variety of such 
algorithms has been proposed and demonstrated based on multidimensional NMR, as 
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described in Z. L. Madi, R. Bruschweiler, R. R. Ernst, J. Chem. Phys, 109, 10603 
(1998). The present invention enables the speeding up of such multidimensional NMR 
quantum computers by several orders of magnitude. This is accomplished by the 
present invention by the method for controlling a multidimensional NMR quantum 
computer with respect to a sample of a material comprising the steps of: (1) rapidly 
monitoring a multidimensional NMR quantum computer with respect to a sample in real 
time by acquiring multidimensional nuclear magnetic resonance spectra within a single 
scan of the sample, and using the spectra for controlling the operation of the 
multidimensional NMR quantum computer. The single scan is carried out with respect to 
the sample and comprises the steps of the methods described herein previously in 
detail. For example, the single scan can be carried out by: (2) partitioning the sample 
into a set of independent subensembles; (3) applying the single scan to the sample by 
exciting the set of subensembles by a time-incremented series of selective excitation 
sequences; (4) generating a signal from each subensemble; (5) homogenously mixing 
the generated signals; and (6) simultaneously acquiring a complete multidimensional 
nuclear magnetic resonance data set for the respective sample. 

Structural elucidations of large molecules, and in particular of systems involving 
proteins and nucleic acids are described in V. Tugarinov, et al., J. Am. Chem, Soc. 124, 
10025 (2002). NMR-based structural elucidations on high molecular weight, complex 
systems will eventually demand the use of a large number of spectral dimensions (over 
4) for achieving sufficient spectral resolution of the peaks. At the same time, however, 
large molecules in general and biological macromolecules under physiological 
conditions in particular, are usually incapable of withstanding the long acquisition times 
hitherto associated with such experiments. By speeding the times required to implement 
such experiments by several orders of magnitude, the present invention bypasses such 
limitation, providing a new way to obtain the structure of macromolecules in their native 
states. This is accomplished present invention by the method for the structural 
elucidations of a large molecule with respect to a sample of a molecule comprising the 
steps of: (1) rapidly monitoring a large molecule with respect to a sample in real time by 
acquiring multidimensional nuclear magnetic resonance spectra within a single scan of 
the sample, and using the spectra for elucidating the large molecule. It may be 
necessary to repeat the single scan a multiple of times to obtain a set of spectra. The 
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single scan is carried out with respect to the sample and comprises the steps of the 
methods described herein previously in detail. For example, the single scan can be 
carried out by: (2) partitioning the sample into a set of independent subensembles; (3) 
applying the single scan to the sample by exciting the set of subensembles by a time- 
incremented series of selective excitation sequences; (4) generating a signal from each 
subensemble; (5) homogenously mixing the generated signals; and (6) simultaneously 
acquiring a complete multidimensional nuclear magnetic resonance data set for the 
respective sample. 

A further subordinate method of the present invention concerns ultrafast 
multidimensional NMR that can be applied regarding in vivo spectroscopy and the 
following of fast metabolic processes, as described in M. S. Cohen, Ed., Physiological 
NMR Spectroscopy: From Isolated Cells to Man, Ann. N. Y. Acad. Scl. (vol. 508, New 
York, 1 987), both in basic research as well as In clinical diagnosis applications, such as 
described in M. A. Thomas, et al., Magn. Reson. Med. 46. 58 (2001). Indeed 
multidimensional NMR spectroscopy on animals and/or humans Is currently hampered 
by the long times that subjects need to reside within the NMR magnet for the completion 
of the experiments, a demand which should be greatly eased by the practice of this 
novel invention. This will enable a new route to the clinical diagnosis of disease as 
described by E. T. Fossel, Diagnosis of cancer by NMR spectroscopy of blood 
lipoproteins and triglycerides. U.S. 4918021 (1990), 18 pp. and by M. A. Brown and R. 
C. Semelka, MRI: Basic Principles and Applications, Wiley-LIss, New York, 1999, 
whether as a pure spectroscopic tool or in combination with methods for spatial 
localization. This Is accomplished by present invention by the method for in vivo 
spectroscopy with respect to a sample comprising the steps of: (1 ) rapidly monitoring 
the sample in real time by acquiring multidimensional nuclear magnetic resonance 
spectra within a single scan of the sample, and using the spectra for evaluating the 
sample. It may be necessary to repeat the single scan a multiple of times to obtain a set 
of spectra. The single scan Is earned out with respect to the sample and comprises the 
steps of the methods described herein previously in detail. For example, the single scan 
can be carried out by: (2) partitioning the sample into a set of Independent 
subensembles; (3) applying the single scan to the sample by exciting the set of 
subensembles by a time-incremented series of selective excitation sequences; (4) 
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generating a signal from each subensemble; (5) homogenously mixing the generated 
signals; and (6) simultaneously acquiring a complete multidimensional nuclear magnetic 
resonance data set for the respective sample. 

Due to the reliance of the invention on magnetic field gradients, further 
subordinate methods concern ultrafast NMR making it also amenable to incorporate a 
variety of MRI protocols. This opens new routes for the accelerated acquisition of clinical 
diagnostic and research MR images, along the lines described in F. Schmitt, M. K. 
Stehling, R. Turner, Echo Planar Imaging: Principles, Technique, Application (Berlin, 
Springer, 1998). The new ultrafast multidimensional MRI invention that results can be 
employed to monitor brain metabolism, pulsating regions (thorax, abdomen), etc. It can 
also aid for the real-time positioning of malignancies and hence as aid in surgical 
procedures. This is accomplished by present invention by the method for magnetic 
resonance imaging with respect to a sample comprising the steps of: (1) rapidly 
monitoring the sample in real time by acquiring multidimensional nuclear magnetic 
resonance spectra within a single scan of the sample at a preselected short time interval 
during the performance of an MRI protocol, and using the spectra for evaluating the 
sample or creating an image of the sample. It may be necessary to repeat the single 
scan a multiple of times to obtain a set of spectra and to obtain a plurality of images. 
The single scan is carried out with respect to the sample and comprises the steps of the 
methods described herein previously in detail. For example, the single scan can be 
carried out by: (2) partitioning the sample into a set of independent subensembles; (3) 
applying the single scan to the sample by exciting the set of subensembles by a time- 
incremented series of selective excitation sequences; (4) generating a signal from each 
subensemble; (5) homogenously mixing the generated signals; and (6) simultaneously 
acquiring a complete multidimensional nuclear magnetic resonance data set for the 
respective sample. 

The extension of single-scan multidimensional spectroscopy can be applied to 
techniques other than NMR. Over the last decade, the concepts underlying 
multidimensional NMR have also been extended and demonstrated in other types of 
spectroscopies: mass spectrometry (A. G. Marshall, F. R. Verdun, Fourier Transforms in 
NMR, Optical, and Mass Spectrometry: A User's Handbook (Elsevier, Amsterdam, 
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1990)), electron paramagnetic resonance (A. Schweiger, G. Jeschke, Principles of 
Pulsed Electron Paramagnetic Resonance (Oxford University Press, Oxford, 2001)), as 
well as a variety of pulsed infrared, optical and UV spectroscopies (S. Mukamel, Ann. 
Rev. Phys. Chem 51 , 691 (2000)). As in the case of NMR these techniques have based 
on monitoring the response of the system as a function of an incremented time variable, 
thereby requiring the collection and processing of numerous individual scans. By 
adapting the ideas of a reversible inhomogeneous frequency broadening coupled to the 
selective excitation of sub-ensembles within the sample, the present invention can 
enable non-MNR multidimensional spectroscopy to be reduced to a single-scan 
acquisition. This is accomplished by present invention by the method for non-MNR 
multidimensional spectroscopy with respect to a system comprising the steps of: (1) 
rapidly monitoring the system in real time by acquiring spectra based on the response of 
the system as a function of an incremented time variable within a single scan of the 
system at a preselected short time interval during the functioning of the system, and 
using the spectra for evaluating the system or creating an image of the system. It may 
be necessary to repeat the single scan a multiple of times to obtain a set of spectra and 
to obtain a plurality of images. The single scan is carried out with respect to the system 
and comprises the steps of the methods described herein previously in detail. For 
example, the single scan can be carried out by: (2) partitioning the system or its region 
of interest into a set of independent subensembles; (3) applying the single scan to the 
system or ROI by exciting the set of subensembles by a time-incremented series of 
selective excitation sequences; (4) generating a signal from each subensemble; (5) 
homogenously mixing the generated signals; and (6) simultaneously acquiring a 
complete multidimensional data set for the system. 

Referring now to Figure 16, there is shown in block diagram a preferred 
embodiment of an NMR/MRI apparatus according to the present invention with the types 
of signals shown in the block diagram shown in the legend box at the lower right of the 
drawing. As shown, the apparatus consists of a magnet 100 that generates a high- 
quality, high-Intensity magnetic field that is uniform within the volume of a sample 106 
undergoing test or monitoring. As shown, the magnet is made of a superconducting 
wire, and includes an ancillary shim coil system for achieving a part-per-milllon 
homogeneity over the desired volume. Included within this embodiment are the coils 



44 



102 for generating the magnetic field gradients required by the invention. The magnetic 
field strength is normally stabilized with the aid of an additional locking circuitry. Within a 
probehead unit 104 in the magnet 100 is positioned the sample 106 and RF coils 108. 

The probehead unit 104 ("probe") contains the sample 106 to be studied, as well 
as ancillary electronic equipment, including the coils 102 needed for generating the field 
gradients required by the invention, the electronics associated with an efficient RF 
irradiation of the spins, and the circuitry for an efficient detection of the spins' signal. 
Numerous such gradient and RF circuits are present in the single probe; the former for 
accounting for the three spatial directions, the latter for the simultaneous Irradiation of 
multiple nuclear species (^H, ^H, ^^C, etc.). 

The gradient wave form generator and driver unit 110 is comprised of a digital 
gradient waveform generator and a gradient driver that translates these digital signals 
into low-level analog currents, which are fed to gradient amplifiers X, Y, Z via lines 140 
where these low-level signals are translated into intense high gradient driving currents 
that are supplied via lines 142 to the gradient coils 102 surrounding the sample 106. 
Three such units are independently present, driving orthogonal x, y and z geometries. 

RF Generation and Irradiation Unit 112 is comprised of an RF unit having a low- 
level synthesizer generating the basic low power RF signal used to irradiate a spin, an 
ampiltude- and phase-control stage capable of creating pulses of different frequencies 
and shapes. Unit 112 is coupled via lines 136 to and feeds high-power amplifiers 114 
that translate these low-level signals into the intense pulses that are fed via lines 138 to 
common lines 139 into the RF coils 108 in the probe 104 for the actual irradiation. 
Common lines 139 serve as multiple RF Channels In/Signals Out. Unit 112 also 
provides a reference RF for the subsequent demodulation of the spins' signal from the 
radio (MHz) to the audio (kHz) range via line 134 to RF Signal Detection/Demodulation 
Unit 124. Several such units are usually present, one per spin species to be irradiated 
during a particular sequence (^H, ^^C, ^^N, ^H, etc.). 
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A Signal Detection Unit receives the RF signal from spins via multiple RF 
channels In/Signals Out line 139 and coupled line 144. The Unit is comprised of an RF 
signal preamplifier/amplifier 120 to effect the requisite preamplification, and then, 
amplification. Preamplifier/amplifier 120 is coupled, In turn, to an RF Signal 
Detection/Demodulation unit 124 including the functions of demodulation and detection, 
which in turn is coupled to an Analog-to-Digital Data Acquisition Unit 126 that contains 
the digitization components, capable of transforming the voltage originally generated by 
the spins following their irradiation into a string of complex numbers (the Free Induction 
Decay or FID). 

A Computer and Display 130 is coupled to Fast Pulse Programmers 132, which 
receive instructions from the computer 130, via bus 150, as indicated in Figure 13. 
These components 130 and 132 are responsible for interfacing to the user, and then 
creating the desired sequence of commands that all remaining units in the apparatus will 
carry out during the course of the operation of the apparatus or the experiment. Thus, 
the computer 130 also provides instructions via bus 152 to the RF generation unit 112 
and to the Gradient Wave Form Generator Driver Unit 1 10. In accordance with standard 
computer practice, user's commands, input via a standard I/O 154, such as a keyboard 
or other such input device, are translated by the computer 1 30 into strings of binary 
digits and logical timing signals, that are in turn executed by the various units. The 
computer 130 containing adequate memory is also usually the final depository of the 
digitized FID, where its data is processed into an NMR spectrum according to the 
algorithms and other information set forth above and display takes place. The Fast 
Pulse Programmers 132 are connected by buses 160 and 162 to provide fast timing 
control to the RF Signal Detection/Demodulation Unit 124, the RF generation unit 112 
and the Gradient Wave Form Generator and Driver Unit 110. 

Referring now to the flow chart of Figure 17, the invention will now be explained 
as it applies to the acquisition of multidimensional NMR spectra. After, a sample 106 
has been placed in the probehead 104 within the magnet 100 in proper juxtaposition to 
the gradient coils 102 and the RF coils 108, the first step S1 is to let NMR magnetization 
generate to obtain relaxation. In step S2, an RF pulse sequence derived from RF 
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generation unit 112 is applied to the sample 106 in combination with a field gradient 
obtained from the Gradient Wave Form Generator and Driver Unit 1 10. In step S3 in this 
fashion a spatially-incremented series of spin states is generated in sample 106, each 
spin state corresponding to different evolution times. This step S3 can be repeated, if 
needed, once per indirect dimension as desired, as indicated by the loop back shown in 
the flow chart from step S3 to step S2. Then in step S4, if needed, a final homogenous 
RF pulse sequence is applied capable of generating an observable signal. 

In step S5 the signals are collected within a single continuous scan, while in the 
presence of oscillating field gradients that reveal the spins original positions (one per 
indirect k dimension), and of a final acquisition time fw. Finally, in step S6, the resulting 
single-scan signal is subjected to a suitable rearrangement along its multiple k axes, 
and to Fourier analysis as a function of the final acquisition time tN. 

The present invention (i.e., system or apparatus described in detail in this 
description of specific embodiments and as generally depicted in Figure 16 or any part 
thereof) may be implemented using hardware, software or a combination thereof and 
may be implemented in one or more computer systems or other processing systems, 
and the capability would be within the skill of one ordinarily skilled in the art of 
programming of computers from the teachings and detailed disclosure provided In the 
foregoing description of the apparatus and the process. The computer system of the 
invention represents any single or multi-processor computer, and in conjunction 
therewith, single-threaded and multi-threaded applications can be used. Unified or 
distributed memory systems can be used. In one example, the system and method of 
the present invention is implemented in a multi-platform (platform independent) 
programming language such as Java, programming language/structured query language 
(PL/SQL), hyper-text mark-up language (HTML), practical extraction report language 
(PERL), Flash programming language, common gateway Interface/structured query 
language (CGI/SQL) or the like and can be Implemented in any programming language 
and browser, developed now or In the future, as would be apparent to a person skilled In 
the relevant art(s) given this description. 
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In another example, the system and method of the present invention, may be 
implemented using a high-level programming language (e.g., C++) and applications 
written for the Microsoft Windows NT or SUN OS environments. It will be apparent to 
persons skilled in the relevant art(s) how to implement the invention in alternative 
embodiments from the teachings herein. 

The Computer system of the invention includes one or more processors and can 
execute software implementing the routines described above, such as shown in Figure 
17. Various software embodiments are described in terms of this exemplary computer 
system. After reading this description, it will become apparent to a person skilled in the 
relevant art how to implement the invention using other computer systems and/or 
computer architectures. 

The Computer system can include a display interface that fonA^ards graphics, text, 
and other data from the communication infrastructure (or from a frame buffer not shown) 
for display on the display unit included as part of the system. 

The Computer system also includes a main memory, preferably random access 
memory (RAM), and can also include a secondary memory. The secondary memory 
can include, for example, a hard disk drive and/or a removable storage drive, 
representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. The 
removable storage drive can read from and/or write to a removable storage unit in a 
well- known manner. 

In alternative embodiments, a secondary memory may include other similar 
means for allowing computer programs or other instructions to be loaded into computer 
system. Such means can include, for example, a removable storage unit and an 
interface. Examples can include a program cartridge and cartridge interface (such as 
that found in video game console devices), a removable memory chip (such as an 
EPROM, or PROM) and associated socket, and other removable storage units and 
interfaces that allow software and data to be transferred from the removable storage 
unit to computer system. 
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The Computer system can also include a communications interface that allows 
software and data to be transferred between computer system and external devices via 
a communications path. Examples of communications interface can include a modem, 
a network interface (such as Ethernet card), a communications port, interfaces 
described above, etc. Software and data transferred via a communications interface are 
in the form of signals that can be electronic, electromagnetic, optical or other signals 
capable of being received by communications interface, via a communications path. 
Note that a communications interface provides a means by which computer system can 
interface to a network such as the Internet. 

The present invention can be implemented using software running (that is, 
executing) in an environment similar to that described above with respect to Figure 14. 
In this document, the term "computer program product" is used to generally refer to 
removable storage unit, a hard disk installed in hard disk drive, or a carrier wave 
carrying software over a communication path (wireless link or cable) to a communication 
interface. A computer useable medium can include magnetic media, optical media, or 
other recordable media, or media that transmits a earner wave or other signal. These 
computer program products are means for providing software to the computer system. 

Computer programs (also called computer control logic) are stored in main 
memory and/or secondary memory. Computer programs can also be received via a 
communications interface. Such computer programs, when executed, enable the 
computer system to perform the features of the present invention as discussed herein. 
In particular, the computer programs, when executed, enable the processor to perform 
features of the present invention. Accordingly, such computer programs represent 
controllers of the computer system. 

The present invention can be implemented as control logic in software, firmware, 
hardware or any combination thereof. In an embodiment where the invention is 
implemented using software, the software may be stored in a computer program product 
and loaded into computer system using a removable storage drive, hard disk drive, or 
interface. Altematively, the computer program product may be downloaded to computer 
system over a communications path. The control logic (software), when executed by 
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the one or more processors, causes the processor(s) to perform functions of the 
invention as described herein. 

In another embodiment, the invention is implemented primarily in firmware and/or 
hardware using, for example, hardware components such as application specific 
integrated circuits (ASICs). Implementation of a hardware state machine so as to 
perform the functions described herein will be apparent to persons skilled in the relevant 
art(s) from the teachings herein. 

Although the present invention has been described in terms of specific preferred 
embodiments, nevertheless modifications and changes will become apparent to those of 
skill in the art from the teachings herein. Such modifications and changes as will be 
readily apparent to those skilled in the art are deemed to fall within the purview of the 
invention as set forth in the claims appended hereto. 
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